Simulations of metastable states near the apex of a force microscope tip interacting 

with an ionic crystalline surface 
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Atoms or pairs of ions picked up by probe tips used in dynamic force microscopy (DFM) can 
be strongly displaced and even hop discontinuously upon approach to the sample surface. The 
energy barriers for some of those hops are of the right order of magnitude to explain the rise in 
energy dissipation commonly observed in DFM measurements at room temperature. The systematic 
computations reported here can explain the infrequent jumps and very low average energy dissipation 
observed low temperature in a previous DFM study on a KBr(OOl) sample. Close to the surface 
we indeed find new states separated by small energy barriers which account for those phenomena. 
These energy barriers strongly depend on details of the atomic arrangement in the vicinity of the 
tip apex. 



I. INTRODUCTION 

Dynamic force microscopy (DFM) has developed into a 
valuable tool not only for surface characterization of non- 
conducting samples, but also for controlled modification 
at the atomic level. This has become possible mainly 
thanks to sensitive measurement modes where the tip is 
oscillated with a constant amplitude in the nm range at 
a resonance frequency of the force sensor [l[. Atomic- 
scale precision is then achieved if the tip apex stays or 
periodically comes to distances at which short-range site- 
selective forces act, thereby causing a measurable fre- 
quency shift. Atom manipulation experiments in that 
mode have inspired computations of changes in the po- 
tential landscape induced by the tip apex and of resulting 
fingerprints in measurable quantities on semiconductor 
surfaces , as well as on ionic crystal surfaces Jf| 0] ■ 
The average dissipation of energy stored in the cantilever 
oscillation also exhibits atomic-scale contrast, even on 
defect-free surfaces, and its magnitude indicates that it 
mainly originates from hysteretic hopping of atoms be- 
tween two or more stable positions @, Q . Sudden but 
infrequent contrast changes, typically more pronounced 
in dissipation images, have been attributed to long-lived 
rearrangements of the tip apex [T^j ■ 

Three different causes of dissipation induced by hys- 
teretic tip-sample interactions must be considered: hop- 
ping on the sample, hopping between tip and surface, 
and hopping on the tip. Hopping on the sample is usu- 
ally not expected because diffusion or other rearrange- 
ments on clean flat terraces of low-index surfaces usu- 
ally involve rather high energy barriers, except for some 
reconstructed surfaces which exhibit bistable configura- 
tions [ri| . The presence of long-lived localized defects 
can be excluded by taking high-resolution images. How- 



ever, in scanning tunneling microscopy studies [12|, as 
well as for DFM on insulating surfaces, mobile adsor- 
bates can merely lead to blurry or streaky images and 
to additional noise in a certain temperature range thus 
causing blurred or streaky images when the scan and hop- 
ping rates roughly match [l3[ . Hopping between tip and 
sample [U, Q , that can in the extreme case even lead to 
atomic chain formation and breaking in some oscillation 
cycles [l4[ , manifests itself indirectly via the average en- 
ergy dissipation. However, in order to unambiguously 
interpret measured results, the third scenario, hopping 
on the tip, must be excluded or else taken into account 
US HI]- This is in general difficult because the struc- 
ture and chemical composition of the tip are unknown. 
For many commonly studied crystals (Si, KBr, NaCl), 
there are indications that sample material is picked up 
by the tip owing to intentional or accidental contact prior 
to or during DFM measurements. Using large-scale sim- 
ulations several groups have characterized force micro- 
scope tips and derived construction principles for realistic 
model tips from comparisons with experiments [1, 0, [HJ ■ 



II. MODEL TIP 

Here, we study possible low-energy configurations of an 
overall neutral KBr tip supporting two additional ions, 
and its interaction with a KBr (001) surface by means of 
extensive computations. The employed code, based on an 
atomistic shell model, was developed for simulations of 
DFM on ionic crystals [ltj E3] and validated in previous 
studies [HI . Such a model tip may represent the end 
of a nominal silicon tip typically used in force microscopy 
experiments decorated by sample material. More pre- 
cisely, our model tip consists of a K + -tcrminatcd cubic 
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cluster of 64 K and Br ions exposing stable {001}-facets 
oriented such that the (111) direction is perpendicular to 
the sample surface. One K + and one Br - ion are added 
near one of the <100> edges meeting at the tip apex, 
as illustrated in Fig. Q] right. The initial configuration 
of these two additional ions is chosen in accordance with 
simulations and a previous experimental study of diffu- 
sion on surfaces of rocksalt-type crystals [2(| ■ This model 
is well-suited to for studying rearrangements of the sim- 
plest moiety likely to be picked upon gentle contact with 
the sample. Alternatively, the supported KBr dimer may 
be the remnant of a broken chain of ions formed during 
tip retraction [3 [2l|. The assumed tip configuration 
is probably more likely than alternative ones involving 
other kinds of defects which produce appreciable force 
hysteresis and energy dissipation in the case of chemically 
similar NaCl model tips [7| ■ The sample was represented 
by a slab containing 6 layers of 10 x 10 ions each. More 
details about the simulation procedure can be found in 
previous publications [H, HU . 

First, the properties of the decorated tip alone were 
studied. One stable and four metastable configurations 
labelled A to E were found by constrained minimization 
while shifting the additional Br~ ion parallel to the edge 
(projected reaction coordinate q) and letting its two or- 
thogonal coordinates and those of all ions in the bottom 
half of the cube relax. The corresponding profiles la- 
belled z = inf arc shown in Fig. [2J In configurations A 
and E, the additional Br - and K + ions are essentially lo- 
cated along the cube edge, and the positions of the Br - 
ion differ between both configurations by approximately 
one bulk lattice constant. The lowest total energy was 
found for the (Br~-terminated) configuration A, because 
the electrostatic field outside the cube is enhanced in the 
vicinity of the low-coordinated edge especially around the 
apex. For the same reason the energy barrier to reach 
configuration A starting from E is much lower than for 
the opposite process. Configurations B and D arise when 
the additional Br~ ion is located above a bridge site on 
one facet adjacent to the cube edge. Between q = 2.1 A 
and q = 2.9 A the initial Br - ion dips into the cubic clus- 
ter while a nearby bromine ion emerges from the cluster 
to form configuration C. Ion exchange processes analo- 
gous to that just described were also found in previous 
simulations of an MgO dimer diffusing on the MgO (001) 
surface |22l|. Similar configurations were found when all 
Br - ions were replaced by K + ions and vice-versa. 

To test the thermal stability of the different tip configu- 
rations, molecular dynamics simulations were performed 
at T = 200, 300 and 500 K for the K+ and the Br~ ter- 
minated tip in configuration A. Below 500 K, no hopping 
was observed over the relatively short simulation time. 
At T = 500 K the former tip showed a transition from 
the E to the A configuration nicely visualized in a movie 
[23l | . In this transition the K + ion did not, however, 
jump directly to the final position. Instead, it moved to 
the position of a nearby K + ion in a configuration similar 
to that called C in the constrained minimizations men- 
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FIG. 1. (color online) Binding energy with respect to an arbi- 
trary offset chosen for visual clarity as a function of decreas- 
ing nominal tip-sample distance for the Br~ terminated tip in 
configurations A and E illustrated on the right. In both cases 
the rigid top half of the tip is approached along the same path 
such that the protruding Br~ in configuration A is initially 
above a Br~ or a K + surface ion. The data calculated above 
K + is offset by -1 eV for clarity. 



tioncd earlier. This exchanged K + ion then moved to 
the position of the K + ion in the A configuration. Such 
exchange processes compete against pivoting around the 
dimer partner as in diffusion on (001) surfaces [2(], [22[ 
and may even be favored in the present lower coordina- 
tion situation. 



III. APPROACH TO THE SURFACE 
CONFIGURATIONS 
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In a second step, the interaction of the tip with the 
sample surface and possible hysteretic processes were 
studied as a function of the nominal tip-sample distance 
defined as the separation of the foremost ions when ignor- 
ing relaxation. Ions in the top half of the tip cluster and 
in two boundary layers of the sample slab were frozen, 
while the rest were allowed to relax. The rigid part of 
the tip was incrementally approached perpendicular to 
the surface such that the protruding Br - in configura- 
tion A was facing a particular surface ion. The results 
were compared to configuration E at the same positions 
of the rigid tip body (Fig. [IJ . If the tip is approached 
above an ion of the same charge, the energy difference 
between configurations A and E becomes smaller, but 
never vanishes. In contrast, if the tip is approached above 
an oppositely charged ion, the energy first drops faster 
then the energy difference decreases and vanishes for the 
K + -tcrminatcd tip, and even becomes negative for the 
Br~-tcrminatcd tip below a critical tip-sample distance 
of about 1 A. One therefore expects the Br~-terminatcd 
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FIG. 2. (color online) a) Energy profiles calculated for the 
initially protruding Br~ tip ion above equally and oppositely 
charged surface ions at a few nominal tip-sample distances. 
The coordinate q of the moved Br~ points along the cube 
edge away from its apex. The data above K + has been offset 
by -1 eV for clarity. All black curves represent constrained 
minimizations at g-increments of 0.1 A. The red points refer 
to subsequent full minimizations, b) and c) Side views on 
the new configuration A*, viewed along the x- and along the 
y-directions, respectively. 



tip to change from configuration A to configuration E be- 
low that distance if the energy barrier between the two 
states can be overcome by thermal fluctuations. 

In order to investigate changes in the energy landscape 
induced by the tip approach, we performed constrained 
minimizations like those discussed in Section |n] at a few 
tip-sample distances. Figure [2] shows that the number 
and character of the metastable configurations changes 
significantly. While at relatively large distances five con- 
figurations analogous to A, B, C, D and E are observed, 
only three remain at close tip-sample distances above the 
Br - surface ion. Lateral hysteresis is also observed, e. 
g. for z = 1 A, indicating the presence of inequivalent 
energy barriers along paths starting from configurations 
A and E. A new configuration (A*), even lower in en- 
ergy than configuration A, appears above the oppositely 
charged K + surface ion for z = 1 A. In configuration 
A* the additional K + and Br~ ions of the molecule are 
bound to both tip and surface. This configuration arises 
when the body of the cubic tip pushes the added dimcr 
aside, until the dimer ions approximately bind to ions of 
opposite charge on the surface as well as on a tip facet, 
as illustrated in another set of movies 24 1 . 



IV. ENERGY BARRIERS FOR 
CONFIGURATION AL CHANGES 

Our study originally focused on hopping between A 
and E configurations, because simulations of diffusion 



on rocksalt type (001) surfaces identified this process 
as the most probable [2(|. The presence of intermedi- 
ate metastable configurations implies that direct hopping 
between A and E is less probable than hopping via the 
intermediate states B, C and D, due to the reduced in- 
tervening energy barriers which are the relevant ones for 
thermal activation. The highest energy barrier between 
any of those configurations represents the bottleneck of 
the process and therefore determines the effective hop- 
ping rate. The bottleneck barrier for the transition from 
E to A ranges from 80 to 175 meV while the bottleneck 
barrier for the opposite transition from A to E ranges 
from 260 to 660 meV for all of the constrained energy 
profiles studied [25[. These values should be compared 
to the thermal energy at room temperature (25 meV) or 
at low temperatures (8 - 40 K 15J, i. e. 0.7—3.4 meV) de- 
pending on the experiment to be considered. Two limit- 
ing cases are of particular interest (l5l . HH : if the hopping 
rate is low enough (one jump every 0.1 — 10 s), individ- 
ual atomic jumps can be experimentally observed. If, in 
contrast, the hopping rate exceeds the cantilever oscilla- 
tion frequency (~ 100 — 200 kHz) the individual states 
involved are averaged over in a dynamic force measure- 
ment, but the energy dissipation due to hops into lower 
energy configurations becomes appreciable. Our results 
imply that, for most potential energy landscapes so far 
considered, the tip configuration would be rapidly flip- 
ping at room temperature on time scales faster than the 
cantilever oscillation, but that tip changes due to hops 
would be frozen out at low temperatures. The main rea- 
son is that in configuration A, the protruding Br - ion is 
subjected to the positive electrostatic potential localized 
around the apex of the cube po| . This lowers its binding 
energy and thus raises the bottleneck energy barrier be- 
tween A and E. Therefore, hopping between A and E can 
only account for infrequent jumps at room temperature, 
but not for the single jumps observed at low temperatures 
0. 

As an alternative model tip one may consider a KBr 
cubic cluster with the additional K + and Br~ ions placed 
in positions of higher coordination such as a facet. How- 
ever, with this premise it is more difficult to realize an 
atomically sharp tip, as required for lattice-resolved im- 
ages without appreciable distortions |27| observed exper- 
imentally. 



V. LOW ENERGY BARRIERS 

The jumps observed at low temperatures might be ex- 
plained by the occurrence of new states similar to A*. 
Indeed, the energy barrier between A and A* being only 
64 meV at a tip-sample distance of 1 A, individual jumps 
could be observed at temperatures between 25 and 30 K. 

In order to investigate the latter possibility in more 
detail, we calculated the total energy as a function of 
tip-sample distance for inequivalent rotations of the Ex- 
terminated tip around the z-axis perpendicular to the 
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TABLE I. Summary of all the energy barriers obtained in 
constrained minimizations above different surface sites. 



sample surface passing through a surface K + site, start- 
ing from alignment of the added KBr dimer with one of 
the <100> axes of the sample at 0°. In this notation 
our previous calculations correspond to 30°. Due to the 
symmetries of the surface and the model tip, full infor- 
mation is contained for polar angles between 0° and 45°. 
At distances larger than about 2.0 A, the data are quali- 
tatively and quantitatively similar for all rotation angles 
(Fig. [3] a) and b)). However, at smaller distances signif- 
icant differences occur. Between 30° and 45°, a jump in 
the energy arises at an angle-dependent critical distance 
below maximum attraction. At smaller distances, the in- 
teraction energy decreases again upon approach. This 
behavior can be identified with the state A* discussed 
above, as evident in movies of the approach J24[. In 
contrast, between 0° and 20°, the energy continuously 
increases further upon approach in the same distance 
regime. The foremost tip ion then remains in a deformed 
state A roughly under the tip apex, while for angles be- 
tween 30° and 45°, the foremost tip atom jumps towards 
a next neighbor surface ion of opposite sign while the 
dimer becomes aligned with a surface <100> direction, 
as well as with a tip facet in configuration A*. For 20° 
the energy versus distance curve is deformed around the 
critical distance, an indication of the proximity to an ad- 
ditional state, but the tip remains in the deformed state 
A. 

We finally studied the stability of configuration A*. 
When the tip is retracted starting from distances larger 
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FIG. 3. (color online) Interaction energy as a function of dis- 
tance for different polar angles of the Br~-terminated tip with 
respect to the [001] axis passing through a surface K + site, 
a) Between and 20° no jump occurs, while between 30 and 
45° (curves offset by -0.3 eV for clarity) hysteresis is observed 
in the investigated distance range, b) Approach-retraction 
hysteresis observed for 30°. c) Results of constrained mini- 
mizations for orientations where no jump was observed and 
(curves offset by -0.3 eV for clarity) for those where hysteresis 
was observed, d) Schematic drawing illustrating two inequiv- 
alent orientations of the tip cluster with the initially on-edge 
adsorbed dimer with respect to a surface unit cell viewed from 
above. 



than critical, the initial values of the energy, force and 
atomic positions are smoothly recovered. Otherwise hys- 
teresis is observed until eventually another jump restores 
the initial energy, force and atomic positions at a larger 
tip-sample critical distance. In Fig. |3]c) an example for 
30° is shown which is also further documented in a movie 
of the simulated retraction [28[. Since our calculations 
are done at zero temperature, the observed hysteresis im- 
plies that the energy barrier between A and A* vanishes 
at the critical distance of approach and that energy is 
gained by jumping to A* at closer tip-sample distances. 
Similarly, upon retraction energy is gained by jumping 
back to state A while the reverse energy barrier van- 
ishes at the critical distance of retraction. Between the 
two critical points a finite energy barrier exists between 
the two states. This variation of the energy landscape 
as a function of tip-sample distance corresponds to the 
scenario proposed by Sasaki and Tsukada Q with the 
modification that the atomic hops themselves need not 
occur in the z direction between the tip and the sample 
but must only be induced by the tip motion perpendicu- 
lar to the sample surface. Results obtained for different 
relative orientations of the tip arc summarized in table 



5 



rotation 


critical 


critical 


hysteresis 


energy 


angle 


distance of 


distance of 


loop 


barrier 




approach 


retraction 


area 


A*^A 


o 


[A] 


[A] 


[meV] 


[meV] 


20 








19 


30 


1.3 


2.1 


200 


91 


35 


1.3 


2.1 


170 


93 


40 


1.1 


2.7 


190 


23 


45 


0.9 


1.9 


90 


40 



TABLE II. Critical distances, hysteresis loop areas and energy 
barriers for different tip orientations. To obtain the data, 
the tip was first approached to the surface up to a distance 
of 1.3 A, then the dimer was forced into state A* by the 
constraint, the system was then fully relaxed and finally the 
tip was retracted. For 20° a configuration similar to A* is 
only reached under the constraint but, upon retraction from 
this configuration, the two added ions remain on the surface, 
so that critical distances and hysteresis are not observed. 



HT1 The energy loss in approach-retraction cycle of the 
tip is equal to the area enclosed between the two distinct 
force-distance curves between both critical distances and 
amounts to up to 0.2 eV which is in the range of what is 
expected from low-temperature experiments jTrij . 

We further characterized the stability of configuration 
A* by studying the energy barrier from state A to state 
A* in constrained minimizations along the previously de- 
fined reaction coordinate q. The tip-sample distance was 
chosen to be 0.13 nm because this is the largest distance 
at which state A* is observed during approach and in a fi- 
nite distance range upon retraction, and so we expected 
that the state could be observed for several tip angles. 
Indeed, between 30° and 40°, where A* was observed in 
energy vs. distance data, state A* is also observed in 
constrained minimizations. For 20°, although state A* is 
not observed during approach, the system can be driven 
into a similar state under the constraint. The energy bar- 
riers range between 19 and 91 meV depending on the tip 
orientation. These energy barriers are of the right order 



of magnitude for explaining the infrequent jumps seen in 
low temperature experiments (l5| . One should, however, 
consider that the chosen reaction coordinate q, was not 
modified to take into account the tip deformation. It 
is therefore possible that even lower energy barriers oc- 
cur in other directions, in particular at close tip-sample 
distances. Obviously, the precise value of the energy bar- 
riers depend on the details of the atomic arrangement of 
the tip. 

VI. CONCLUSIONS 

We conclude that ion pairs picked by force microscope 
tips decorated by sample material can be strongly dis- 
placed upon approach to the sample surface, in particular 
at close tip-sample distances where they become bound 
to both tip and sample. The resulting hops have com- 
ponents parallel to the sample surface. Some of those 
hops can account for rapid nipping of the tip configura- 
tion at room temperature. In addition, some metastable 
states which occur at close tip-sample distances can be 
separated by energy barriers that are low enough to ex- 
plain infrequent individual jumps observed at low tem- 
peratures. The adsorbed dimers have a tendency to align 
with ions of opposite charge on the sample surface. Ion 
exchange processes previously identified in a study of sur- 
face diffusion are preferred for some of the investigated 
hops. Back and forth hops between metastable configu- 
rations result in a hystcretic force as a function of dis- 
tance. The energy dissipated by such hops is in the range 
of what is expected from low-temperature experiments 
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